SNUTP 01-023 



Late decaying axino as CDM and its lifetime bound 

Hang Bae Kim'** and Jihn E. Kim^'f 

"'Department of Physics, Lancaster University, Lancaster LAI 4^B, UK 
^Department of Physics and Center for Theoretical Physics, Seoul National University, Seoul 

151-747, Korea 

Abstract 

The axino with mass in the GeV region can be cold dark matter(CDM) in 

the galactic halo. However, if R-parity is broken, for example by the bilinear 

terms /ia, then axino(d) can decay to v + ^. In this case, the most stringent 

bound on the axino lifetime comes from the diffuse photon background and 

we obtain that the axino lifetime should be greater than 3.9 x lO^^f^a/i s 

which amounts to a very small bilinear R-parity violation, i.e. /Xq < 1 keV. 

This invalidates the atmospheric neutrino mass generation through bilinear 

R-violating terms within the context of axino CDM. 
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The observed rotation curve of the halo stars |jl|] requires to fill the galactic halo with cold 
dark matter(CDM) such as axion 0-0, the lightest supersymmetric particle(LSP) the 
axino LSP [§,0, and wimplzilla ^j. The axion is motivated from the solution of the strong 
CP problem a la Peccei and Quinn(PQ) |0, which is required to be very light |]T0|. The 



LSP is motivated from R-parity conservation in the supersymmetric solution of the gauge 
hierarchy problem, where R-parity is defined as (_i)3B+i-25_ p^j, ^^le LSP to be CDM, its 
mass is around 100 GeV The wimpzilla is 10^^^^^ GeV stable particle. 



The possibility of axino dark matter, which is our interest in this paper, has been sug- 



gested from time to time as the hot DM, the warm DM and the cold DM possibilities |]12 
Theoretically, it arises in SUSY theories with a spontaneously broken PQ symmetry. In 
supergravity it is necessary to have a period of infiation to dilute the very weakly interact- 
ing gravitinos. But it is thermally produced in significant numbers even after the infiation, 
which requires a low reheating temperature, < 10^ GeV |jl3|. A similar study for axino 



requires a much lower reheating temperature of order TeV in which case O(GeV) axino 
mass is allowed. Then axino can be a DM candidate Here, we focus our attention on 

this CDM axino. Since there is no reliable constraint on the axino mass [0,0, the axino 



is assumed to be the LSP. Then, the most important question is whether the LSP (axino) 
is absolutely stable (due to R-parity conservation) or unstable. Since the stable axino case 
has been extensively studied [0, we restrict our attention on the R violating case. 

If a CDM candidate is proposed, it is of utmost importance to devise a scheme to prove 



its existence experimentally as in the cases of the very light axion ||T5| , p!6| and the LSP 0,|TT|. 
The axino CDM lacks this kind of possible detection mechanism due to its extremely weak 
interaction strength if R-parity is conserved. However, if R-parity is broken, it may be 
possible to detect its decay products. Therefore, it is worthwhile to study possible decay 
mechanisms of CDM axino. The detection possibility by the axino decay relies on the axino 
lifetime around > 10^^ s after the galaxy formation era. 

The R-parity conservation seems to be an attractive proposal for proton stability. How- 
ever, R-parity is not dictated from any deep theoretical principle. For example, if there 



exists an SU{3) x SU{2) x U{1) singlet superfield which is needed from the see-saw mech- 
anism LaNH2 where is the lepton doublet of the a-th family and Hi is the Higgs doublet 
{i = 1,2), we may write a renormalizable superpotential, Wn = mN"^ + fN^. From the 
coupling to the observable sector fields for the Dirac mass coupling, is required to carry 
-1 unit of the R-parity quantum number. But this R-parity is broken by the A^^ term. Also, 
A^ obtains a vacuum expectation value — 2m/3/ which is too large for neutrino phenomenol- 
ogy. Namely, R-parity conservation is not guaranteed a priori at the SM level. Thus, if a 
singlet A^ is introduced, one should impose (at least an approximate) R-parity conservation, 
namely we should impose / = in this example. If R-parity is broken, it must be done so 
very weakly. In this paper, for simplicity of the discussion, we restrict our attention on the 
bilinear R-violating terms, 

fiaLaH2 (1) 

where a = 1,2, 3. fi^ is bounded by the eV order neutrino mass. Without an explicit 
statement, this bound applies to the heaviest SM neutrino, presumably the tau neutrino. 
With the R-parity violation, the u^- mass arises from the see-saw type diagram with an 
intermediate zino line with two insertions of the R-parity violating (u). Also, it can get a 
contribution from the intermediate H2 line with two insertions of /is. These give similar 
conclusions and we discuss the fj,^ case for an explicit illustration. Then, /X3 is bounded as 

|/i3| < M'Jl^^^ MeV (2) 

where M^o rp^y is neutral Higgs mass in units of TeV. In this paper, we introduce dimen- 
sionless numbers: for a small coupling e, e_„ represents it in units of 10~", for MeV order 
masses niMeV in units of MeV, for GeV order masses mcev in units of GeV, for large mass 
m\n\ in units of n GeV, and for super large mass F12 in units of 10^^ GeV. 

The bilinear R-violating parameters have been extensively discussed in regards to the 
neutrino oscillation |]T^ within the above bound (^. In this paper, we will draw a conclusion 



that this is not consistent with the CDM axino, already from the observed diffuse gamma 
ray background. 



With the bihnear R-parity violation, we expect the following decay modes of the axino 

a —> 1/ + 7(or /"'"/"), a — i> z/ + a, a— >r"'" + 7r~, etc. (3) 

where rria > rrir + is assumed. 

To estimate the partial decay widths, let us assume the following R violating axino 
interaction 

C-a-dec^y = cq ^o?/' , or iei^^F^''d'j5['j^, 'Ju]i^ (4) 

where Fa is the axion decay constant (including the division by the domain wall number 
Now), Ffj^u is the field strength of a spin-1 field A^, is a scalar field and ip is a. fermion 
field (Dirac or Majorana field). Then, the lifetime of axino becomes 



iev ■ i^-^K-n ■ Po, or 2.57 x lO-^e^F,;,^^ ■ P,)'' [sec] (5) 



where = 1, 2, respectively for the Dirac and Majorana ip, we neglected m^, and Po,i are 
phase space factors. For a massive final fermion, PqTi ~ ^ '^V'/"^^)~"'^(l + f^ip / f^a)''^ ■ Let 
us proceed to discuss several possibilities of O(GeV) axino decay. 

Firstly, for the a ^ v + a decay, we note that the axion multiplet couples to the standard 
model chiral fields, below the PQ symmetry breaking scale, as exp(2(5A/-Fa)W^Q where Wq 
carries —Q units of the Peccei-Quinn(PQ) charge, and A is the axion supermultiplet. The 
heavy quark axion models do not allow these tree level couplings since the SM fields are 
neutral under the PQ symmetry, but the lepton-coupling type models can lead to this 
kind of couplings. Our interest is on the superpotential for the axino-neutrino coupling, 
myVi' e-x.Y>{2iQjjA/ Fa) where Qi, is the PQ charge of the neutrino, and ^aLaH2e^^^f^'' . Set- 
ting Qy = 1/2, Q = 1/2, we obtain W = rriiyUue^'^^^'^ and HaLaH2e^'^^^'' from which we 
obtain the relevant terms for the axino decay —{m^A'^ /2Fa)i'i' and —{^aA'^ /2Fa)LaH2, re- 
spectively. The R-parity violation by the vacuum expectation value of a sneutrino, va = (u), 
or by allow the following Yukawa couplings 

^d^ua = I or -j^ I X aua (6) 



from which we estimate eo — 10~^^mu^evVu,Gev/ F^i2 fJ'a,MevV2,[ioo]/ u, respec- 

tively, and V2 = {H2). Thus, in view of Eq. @ eg is too small (i.e. ~ 10^^ s) and the 
decay mode a u + a is not important cosmologically. 

Second, the decay d r+7r+ occurs through the bihnear R-parity violating term fiaLaH2 
with a = 3, given in Eq. (|I]), which allows the mixing of f and H^. Then, the coupling 
{mT-/Fa)f^Td and the Yukawa coupling ~ {nid/vijqd'^Hi give an effective interaction @ 
with 



_ ni^Uml 1 1 



where AM^py is the f-H^ mixing parameter. The effective interaction (|^) with ip = t and 
(j) = TT^ arises from the tree diagram with the f — H~ mixing insertion in the intermediate 
scalar propagator with the four external fermions, a, r, d, and u. In estimating eo, we used the 
PCAC relation in obtaining the matrix element (OldRULln^) ~ fT^m^/md. The superpartner 
masses for the gauge hierarchy solution are around 100 GeV. The bound on the stau-charged 
Higgs mixing parameter AM^py is bounded from the tau neutrino mass bound, mi, < 1 
eV. For the R-parity violating bilinear coupling fi3L^H2, the mixing parameter is estimated 
as AM|py = 2fi*^fi. Using Eq. (D, we obtain AM|py < GeV^ Then, we 

estimate cq < 3.71 x 10~^^[cos/3]~^F~]'^2/^TeV'^j^- [ioo]^^[^ioo]"^fitf the rvr decay mode, 

and the axino lifetime must satisfy a bound 

r, > 0.959 X 10- [sec] ■ P^"^ (8) 

2,1 eV 

where the MSSM parameter tan/? = V2/V1, and we assume that Pq is nonzero, i.e. the 
O(GeV) axino has mass > 1.92 GeV. 

Third, we note that the interaction vrfd, arising from stau intermediate state and R- 
parity violating insertion of neutrino-!?/ mixing, is not important. 

Finally, we note that the decay a — >^ z/ -|- 7 (or occurs through the anomaly term 

i, 

= ■ ^ ■ ^.75[7^ r]aF,. (9) 

Ibvri'a fi 
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where Ca-y-y is the axion-photon-photon couphng which depends on models [|T8[, and F^y 
is the photon field strength, and the photino-neutrino mixing parameter c'/iaZ/i has been 
introduced. Note that c' is 0(<1) and yU^ is 0(<MeV). It turns out that this photon mode 
constitutes the most important contribution in the axino decay. The gluon anomaly term 
can be considered, but it is not important since we must consider intermediate gluino and 
squark lines. From the interaction Eq. (H), we estimate 

ei = 1.99 X 10"'^Ca77c'/ia,A4'ey/i[ioop (10) 

giving the axino lifetime 

Ta - 9.8 X 10^^ [sec] C^^V"Va,MeV/^[100]^a,12"^S,Ly (H) 

The RHS of Eq. (Il^) can fall in the cosmologically interesting scale for Fa slightly smaller 
than 10^2 GeV and = O(GeV). 

This leads us to the estimation of the cosmological abundance of axino. The decoupling 
temperature of axino is of order the PQ symmetry breaking scale p!^,|7|]. Thus, for an O(GeV) 



axino, inflation must end below the PQ symmetry breaking scale so that axinos produced 
at the decoupling temperature is sufficiently diluted. However, if the reheating temperature 
after infiation were high enough, a significant number of axinos would have been reproduced 
thermally and can constitute cold dark matter. Here, we are interested in this thermally 
produced axinos after infiation. In this scenario, the number density depends on axino mass 
and reheating temperature Tr. For O(GeV) axino in R-conserving theories the reheating 
temperature bound is Tr < 100 — 1000 GeV from the condition that the thermally produced 
axinos do not exceed the critical energy density as estimated in Ref . . On the other hand, 
if there exists R-violating terms, then the lightest neutralino decays to l^l^v within 1 s 
for /ia/yU < 10~^, which occurs through the diagram x\ ~^ W^X^ ^ind the mixing of 
and l^. Namely, a neutralino predominantly decays to SM particles, which is harmless at 
a later epoch. Thus, to produce O(GeV) axino copiously in R-violating theories, reheating 
temperature must be raised. 



In the present estimation of reheating temperature Tji, we ignore the gluino decay to 
axino since R-parity is broken. Taking into account the photon thermahzation after the 
decay of MSSM SUSY particles and e+e^ annihilation, the axino energy density at present 
is PaiT^) = 'n^a(27r^/45) ■ (43/11) ■ T^Y^^ {Tr) where the thermal production of axino Y^^ 
is calculated just from the scattering processes. Thus, representing = x (the critical 
energy density), we have maY^^ ~ 0.72 eV(fia/iV0.2). For Fa ~ lO^^ GeV, Via ~ 0.3 and 
O(GeV) axino, Y^^ ~ 0.5 x 10~^ and we can read Tr at around 200 GeV from Fig. 1 of 
Ref . . In this region, if one considers the gluino decay in R-conserving theories, Tr should 
be a factor ~ 2 smaller. The neutralino decay in R-conserving case is not important in this 
region, but can be very important for Tr < 100 GeV 0. 

The present dark matter density in the galactic halo requires a significant amount of 
dark matter. In the present case, there is the other candidate for dark matter, the axion. 
For the axion CDM, we are restricted to Fa ~ 10^^ GeV. But with the CDM axino. Fa can 
be lower as far as Fa > 10^ GeV. If the axino lifetime is of order the age of the universe, 
then there remains a significant number of axinos which can be detected in experiments in 
search of proton decay. The interesting decay mode for axino detection is the mode. 

For 10^ sec < < tree where tree is the time at the recombination, there can be an 
allowed region Ta > tmin so that the decay products are not copious enough to dissociate 
the light nuclei. After the time of recombination, photons are not effective to scatter off the 
neutral particles, and hence tmin < ^rec- Since we are considering the axino lifetime > 10^^ 
s, the late decaying axino is safe from destroying light nuclei. 

If axino decays to z/ + 7, the underground neutrino detectors can detect the photon. 
The photon energy(= Fry) of 10 GeV, i.e. axino mass of 20 GeV, is the boundary for using 
different search types for the Cherenkov rings. If < 10 GeV, the Compton scattering 
on an atomic electron kicks out a high energy electron whose Cherenkov radiation can 
be detected. If E-y > 10 GeV, then e'^e~ pair production off nucleus dominates and the 
Cherenkov rings from these pair can be detected. 

From the super-K detector, one establishes the proton lifetime bound of 10^'^ seconds 
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19| . These detectors use baryons in water, hb = Na/ctci^ where A^^ = 6.023 x 10^^. On the 
other hand axino as CDM now has the local number density of order = O.Sm^Q^y/cm^, 
giving the ratio Ua/riB ~ 5 x 10"^^m^Qgy. If we require the detection rate of axino decay 
the same as that of proton decay, we obtain [Tpna/TaUp) ^ 1. Thus, axinos are detectable at 



the rate of counting proton decay debris with proton lifetime of 10^^ seconds [|T^ if 

1.6 X lO^*^ , , 
r-a ^ [sec]. (12) 

fna,GeV 

Since we have not observed this kind of events, we obtain > 1.6 x lO^^m^Q^y s. 

However, the most stringent bound comes from the diffuse gamma ray background. The 
classical study on this effect has been published more than 20 years ago 0]. The observed 
flux JF^ is bounded by ^1 



^ < (10^3 ~ 10-^) Ealv cm-^sr-^sec-^ (13) 

where E is the decay photon energy at present. The figure 10"'^ is for the conservative 
bound applicable to the whole observed range of E. For E = IMeV ~ lOGeV where we are 
interested, 10"^ gives a good fit to the data. 

On the other hand, the decay of axinos produces the diffuse photon flux at present, for 
0< E < ma/2 



where Eq = (ma/2)(ra/to)^^^ and to is the age of the universe (~ 4 x lO"*^"^ s). 

For tj-ec < Ta < ^0) niost axiuos have decayed and the flux has a peak at E = Eq with 
a value ria/STce. Assuming the critical axino density, we have a condition that the photons 
from axinos decay do not exceed the observed flux, 

(Ta/to)'/' < 1.4 X lO-'Qah^ (15) 

which is inconsistent with the condition > tree (~ 10^^ s). 

For Ta > to, axino decays are increasing at present and the maximum flux is 
(3na/87r)(to/Ta)e~*''/^'^ at E = 'ma/2. Comparing this with Eq. (|T3l), we obtain 
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Ta,sec > 3.9 X lO^^-'^fia/l. (16) 

For the region satisfied by Eq. (|16]), the bound on /Iq, is very stringent, i.e. less than 
O(10^~^ eV), hence the idea for neutrino mass generation via bihnear R-parity violation is 
not consistent with CDM axino. A very conservative upper bound on /Iq, is 1 keV. 

In conclusion, we searched for the detection possibility of axinos as CDM with Tr ~ 200 
GeV. The diffuse gamma ray background gives a very strong bound on bilinear R-parity 
violating parameter Even if /Iq, is of order keV, it can be detected by diffuse gamma ray 
background observation. On the other hand, with O(keV) fia SUSY generation of neutrino 
oscillation parameters through bilinear R-parity violation is not achievable with CDM axino. 
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